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INTRODUCTIOK 


The general systematics of nuclear binding energies are well parameterised 
by the semi empirical mass formula. The volume term in this formula 
predicts that, in absence of coulomb interactions, an infinite, homogeneous 
fluid of equal numbers of neutrons and protons will be self bound with a 
binding energy per nucleon of Eg © 1620.5 MeV. The symmetry term in the 
semi empirical mass formula predicts an energy cost of 32 MeV/nucleon in 
moving away from neutron-proton symmetry. Measurements of the central 
densities of heavy nuclei suggest that the infinite nuclear matter fluid 
will saturate at a density p, © 0.17+0.01 nucleons per fm* corresponding to 
a Fermi wave number kp * 1.3720.03 fm’. While it is difficult to 
unambiguously separate bulk and surface effects, the observation 
"breathing mode' excitations at 12-15 MeV excitation energy in nuclei 


suggests a bulk compressibility modulus for the nuclear matter fluid 
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k = 997 ŻĘ | a 250 + 50 MeV (1) 
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Attempts by nuclear many-body theorists tu explain this data have a long 
and frustrating history which we shall review in this short communication. 
BRUECKNER AND THE COESTER LINE 
The simplest model for nuclear matter is a fluid of non relativistic 


neutrons and protons interacting through a two-body potential. The 


hamiltonian is then 


Hap aaa 2 p owyj (2) 
i 2m i>j 


The two-body potential must fit the nucleon-nucleon scattering phase shifts 
and the bound state deuteron properties. There are a number of such 'phase 
equivalent’ potentials on the market. They are either completely 
phenomenological in nature or guided to a greater or lesser extent by meson 
exchange models. In general these potentials have a long range 
one-pion-exchange tail, an intermediate range at which they are attractive, 
and which is frequently attributed to an enhanced two-pion-exchange, and a 
short range repulsive core, which may be attributed to vector p- and w- 
meson exchange. There is clear evidence for the non central nature of the 
interaction and the observed quadrupole moment of the deuteron would appear 
to require a tensor force. although there is continuing debate as to the 


strength of the tensor term. 


whatever the details of the interaction, the strength of the repulsive core 
ałone precludes a direct application of the Hartree-Fock prescription for 
handling many fermion systems. In the Hartree-Fock prescription a trial 
wave function in the form of a Slater determinant 


4 = (A!)7% deto;(x,5) (3) 


is chosen and the expectation value of the hamiltonian (2) is optimised 


with respect to the single-particle basis ©;. This leads to the equations 
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for the single-particle basis, and for the ground-state energy expectation 


value 


? 
Eo $4 [eit [ey Jory ag (5) 


The tria] wave function (3) describes a system of particles whose motion is 
correlated with that of its neighbours only through a mean field. In the 
nuclear matter problem the strong short-range forces induce two-body 


correlations which cannot be ignored. 


Keith Brueckner attempted to account for those correlations by allowing 
pairs of particles to undergo multiple scatterings and summing the so 


called ‘ladder diagrams', thus generating a reaction matrix. 


<ij]G|1j> = <ijjV|ij> „ J<ijfvjab>]* < e 
5G} 43 jivlij> + ab eat 


The prime on the summation in eqtn (6) denotes that the intermediate states 
a, Op must honour the Pauli principle and be unoccupied in the 
uncorrelated ground state © of eqtn (3). The energy denominators of eqtn 


(6) can be written 


Eabij = at €b ĉj - èj sy 


The basis of single-particle states occupied in o is now given by 


2 
- Fa 71 908.) + g f oy" (aed GQ jka) Foy (asdoglaa) - og (80g (3a) da 
"MS oj (x,) (8) 


and the lowest order expression for the ground-state energy expectation 


value is 


2 
EU) = x r fes" w [- = w+ ej] oj(x)dz (9) 
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Because of the structural similarity of eqtns (8) and (9) to eqtns (4) and 


(5). this approach is called the Rrueckner-Hartree-Fock prescription (BHF). 
It must be remembered however, that the true Hartree-Fock prescription is 
bused upon the variatjona) principle while the BHF prescription simply 
represents the lowest order term in a partially summed perturbation theory 
expansion. The discerning reader wi]] have noted that I have not fully 
defined the BHF approach because 1 have not defined the basis of 
single-particle states unpopulated in % which are required in the summation 
of eqtn (6). The only safe choice is to use a system of plane waves 
orthogonalised to the states occupied in % and to calculate the single 
particle energies ea and ep as expectation values of the kinetic energy 
operator. A justification of this approach may be that the interaction V 
is so strong that it scatters particles so far out of the Fermi sea that 
their potentia] energy is negligible compared with their kinetic energy. 


Going beyond lowest order perturbation theory it is conventional to 
rearrange the series so that it is no Jonger in order of ascending powers 
of the reaction matrix but rather in ascending order of the number of 
vacancies in the Fermi sea - the so called 'hole-line' expansion. For a 
central, repulsive hard-core potential of radius rę the hole-line expansion 


can be shown to converge with powers of the small parameter 


x = (Tc/rg)* (10) 


where the density of the nuclear fluid is 


p za L/int) O (11) 


At the nucledr saturation density po we have Tc/ro = '/, and converge of 


the hole-line expansion would seem assured. 


The results of BHF calculations for a wide range of phase equivalent 


potentials are summarised in Figure 1. 


The results for a wide range of phase equivalent potentials scatter about 
the line AB, known as the Coester line, and most certainly do not account 
for the semi enpirica] mass formula data represented by the black triangle. 
It is possible to choose a nucleon-nucleon potential which will saturate at 
the correct density but this fails to give sufficient binding energy. It 
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is possible to choose a potential which will provide sufficient binding 


energy. but then it wi]) saturate at too high a density. 
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Fig 1 The saturation of nuclear matter. The black rectangle represents 
the predictions of the semi empirical] mass formula. The line AB 
is known as the Coester line. The results of BHF calculations for 
nuclear matter with a wide range of phase equivalent potentials 
scatter about this line. 


A closer examination of the results of BHF calculations shows that it is 
the tensor force which plays an important role in producing saturation. 
Thus it is no surprise to learn that those bare nucleon-nucleon 
interactions which have a strong tensor force saturate first and yield 
results at the upper end A of the Coester line while the results with 


interactions containing weak tensor forces cluster towards the lower end B. 
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JASTROW CORRELATIONS: BRUECKNER JUST: PIED? 


Given the complexity of a fully self-consistent BHF calculation, 
particularly if there is need to go beyond lowest order as might be 
required at supernuclear densities, eg. in order to describe the core of a 
neutron star. it is attractive to explore alternative calculational 


procedures. 


One possibility would be to stay with the variational principle and improve 
the choice of tria] wave functions. Instead of the Slater determinant (3) 


consider 
P = FO (12) 
where F is a symmetric product of two-body correlation functions 


Fe s f (13) 
1>j ij 


designed to accommodate the effect of particularly strong components of our 
nucleon-nucleon interaction. In the Jastrow approximation -the fij are 
central functions of the separation between the pair of particles which 
reduce the two-body wave function at short distances, where the repulsive 
core forces the particles apart, and which heal to unity at large distances 


where the interactions are extremely weak. 


The full wave function (12) is extremely difficult to work with and it is 
usua] to make some form of cluster expansion of F based upon the relatively 
low density of the nuclear matter problem as measured by the small 
parameter «x (10). The probJem “then is that, if the cluster expansion is 
terminated at any finite order and the correlation functions fij freely 


varied, the energy functional 
E([F])= <*]H |©>Ąpprox (14) 


will not saturate, the so called Emery problem. The notation < >approx 
means that the expectation value in some limited cluster expansion 
approximation. This problem can be traced to the fact that, In searching 
for a lower bound to the energy in the variational calculation, one can 


increase the radjus of the wound in the correlated two-body wave function 
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to completely eliminate the repulsive core and at the same time reduce the 
rate of healing of the correlation function so that the kinetic energy is 
negligible. Then we can adjust the long range part of the correlation 
function to give unbounded weight to the attractive tail of the 
interaction. Of course, such long range correlations are inconsistent with 


a termination of the cluster expansion. 


What is required in the variational calculations is a constraint on the 
correlation functions fij which wil) ensure the convergence of the cluster 
expansion. Then for a constraint 


g ([f]) = 0 (15) 


the variational functional is 


s | | <elt lm approx - X8 (igh) } = 0 (16) 


with à an undetermined Lagrange multiplier. A natural choice for the 
constraint is that it should retain the convergence properties of the 
hole-line expansion. This can be achieved by requiring that on average the 
wound volume around any given nucleon contain at most one further nucleon, 
thus making the restriction to two-body clusters very plausible. Such 
calculations are generally called lowest order constrained variational 
calculations (LOCV). 


When LOCV calculations are performed for a wide range of phase equivalent 
nucleon-nucleon interactions the results agree to an impressive extent with 
BHF calculations using the same potentials. The Coester line is reproduced 
and once more the saturation predictions disagree with the semi-empirical 
mass formula. 


We note that both LOCV and BHF calculations not only fail to reproduce 
nuclear matter saturation but that they predict very high 
compressibilities, typically K * 400 MeV which is inconsistent with the 
observation of breathing mode excitations at 12-15 MeV in nuclei. 
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TENSOR CORRELATIONS: THE MYSTERY, AND THE BINDING-ENERGY, DEFPENS 


The full variational] principle states that the results of some restricted 
varialiona] calculation wil] be an upper bound to the true ground state 
energy and that the calculated ground state should approach the hamiltonian 
ejgenstate as the restrictions on the scope of the variationa] calculation 
are removed. Within the class of LOCV calculations the simplest relaxation 


would be to free the form of the two-body correlation functions fij: 
The form of our two-body potentials is 


A A 
Vij r v ij (17) 


where in different reaction channels » we have the freedom to include 
central, tensor, spin-orbit and quadratic spin-orbit forms as appropriate. 
It seems natura) to accord the correlation functions the same deprees of 


freedom and write 


fij = Ę fij 7 (18) 


with 
fig. = fe>(rgj) + fq>(rgz) Sig + fus*(rig) LS + foòlrij (ES)? (19) 


For all singlet states there are no tensor or spin-orbit terms and thus à 
is simply the partial wave L = 41 - Sy, for the triplet states with J=L 
again à is simply the partial wave but triplet state with J = L+1 are 
coupled to the partial wave L' = L+2 states by the tensor force, eg, 
15,-7D,, °P,-?F, etc, and in these cases X refers to the two orthogonal 


reaction channels associated with the pairs of partial waves. 


Since the new correlation functions are designed not simpły to take into 
account the wound in the two-body wave functions produced by the repulsive 
core ‘of the two-body interaction but also any other possible strong 
correlations induced by other elements of the force, the constraint used in 


the Jastrow cluster expansion no longer seems appropriate. 


The Nuclear Matter Saturation Problem 257 


n a M o 
The single-parlicle energies e; in the original Natree-Tock approach (4) 

are themselves Lagrange multipliers chosen to assure the orthonormality of 
the single particle functions oy and hence the normalisations of the tria] 
many-body function % (3). In this we shal] use as our constraint (15) the 


normalisation of our two-body correlated wave functions. 
Wj = fqj (2!)7% det o; (%5). (20) 


Later we sha)) check the convergence of our cluster expansion and show that 
our termination in lowest order is valid. One consequence of our change in 
constraint is that because jj contains a wound at short distances, then, 
in order to be normalised, fiy must overshoot unity in a region where the 


interaction is attractive and hence we shall get a slightly increased 


binding energy. 


It has been known for many years that in perturbation theory the tensor 
force converges relatively slowly. This was recognised by Kuo and Brown 


in their analysis of the effective central *S, interaction in nuclei where 
they wrote 


E (21) 


g 3 > z 
<*s,|veff|35,> = <°S, [Vg] ?S.> + z <"*S IV = z n=D,> i 


Using an approximate closure argument 4*S,|S,ż|?5,> = 8 they obtained the 
result 


eff 35 
Vag $ Ve " + BVq*/Eętę (22) 


Kuo and Brown reasoned that Epgęg was 200-300 MeV and that this second order 
tensor contribution yielded a considerable increase in the attractiveness 
of the °S, effective interaction. By putting tensor components directly 
into our trial wavefunction we include these effects and many more. It 
transpires that it ts the tensor components in the correlation functions 
(19) which are the really significant new feature and we shall not refer 


further to the spin-orbit components. 


When the new LOCV calculations are carried out for the same range of phase 
equivalent potentials the results cluster in the cross hatched box in 


figure 1. Now the results for various phase equivalent potentials are very 
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similar, instead of being spread along the Coester line they are al] 
grouped together. This is not surprising, we would expect the effect of 
tensor correlatious to be largest for those interactions with strong tensor 
forces, je those which were at the end A of the Coester Jine, and least 


for those with the weakest tensor forces. 


Recently perturbation calcuJations based on the hole-line expansion have 
been extended to fourth and fifth order confirming the results of these 
generalised IOCV calculations. The bad news is that the results of the 
cross hatched area do not overlap the empirical black rectangle. In 
addition the compression modulus remains uncomfortably high with 


K > 350 MeV for most of the phase equivalent potentials. 


Perhaps the most alarming aspect of the results is that the variational 
calculations predict binding which is greater than the empirical] value. 
This phenomena has happened at least once before in the history of physics. 
As variational calculations for atomic structure calculations were refined 
there came a time when the Schrodinger equation for the helium atom yielded 
a ground state energy below that which was experimentally observed. This 
violates the upper bound condition of the variational principle. The 
solution to this apparent paradox is that the model being used must 
incompletely describe the physics being investigated. In the case of the 
helium atom the breakdown was traced to the inadequacy of the non 
relativistic Schrodinger equation. In-the next section we shal] explore 


the missing element in the nuclear matter problen. 


ISOBARS TO THE RESCUE 
Let us review our mode] which consists of a fluid of non relativistic 
nucleons interacting through pair potentials. The potentials necessarily 
contain tensor forces and these we know can be generated by the exchange of 
pions. So at root our model is a soup of nucleons and pions. The 
w-nucleon scattering cross section is dominated by a broad resonance at 
~“ 300 MeV. This fs variously referred to as the 3-3 resonance, the 4 
particle or the N*(1232) isobar. Thus it would seem natura] to generalise 
our trial wave function to allow for the possibility that our nucleons may 
become excited into N* isobar states. This economically achieved by 


extending our correlation function (19) 
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fi Ty © E EE T S er N s (23) 
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where Ean is the generalised tensor operator connecting the two-component 
nucleon spinner to the four component spin ?/, isobar. Such a term arises 
first in the coupling of the 'S,(NN) channel to the DUNN”) channel. The 
analogy with the tensor correlations is now complete and is exhibited 


graphically in figure 2. 


Fig 2 a) The intermediate D-state admixture in the two nucleon °S, 
channel induced by a tensor force, eg, one-pion-exchange. 

b) The intermediate D-state N-N* (represented by the heavy line) 

induced by a one-pion-exchange dugeract ie in the *S, 


nucleon-nucleon channel. 


Naturally our hamiltonian must now be extended to include the isobar 
channels and this is done using the N-N* potentials of Green, Sanio and 


Niskanen based on the w-N production of N* as analysed by Arenhove]. 


The contribution of figure 26 to the N-N 'S, scattering amplitude is 
subsumed into the phenomenological potential fits and is responsible for a 
significant part of the attraction in this channel. Note that 


one-boson-exchange models of the N-N interaction which try to mimic the 
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two-pion-exchange with a single fictitious o-meson exchange miss the 
richness of this phenomena and thus dangerousły hide real physics. When 
two nucleons interact in the nucleus the amplitude in figure 2h must be 
modified, the range of momenta allowed to the intermediate nucleon being 
restricted by the Pauli principle to Jie outside the Fermi sea. Thus we 
must subtract the contribution of this amplitude for all intermediate 
nucleon momenta k<kp. This subtraction turns this attractive amplitude 
into a repulsion. Such an effect had been investigated many years earlier 
by Green and co workers who found it to be a small effect or order 1-2 MeV 
at the empirical] saturation density. Starting from point A on the Coester 
line this further increased the discrepancy between BHF calculations and 
the seni empirica] mass formula. Now we see that the repulsion increases 


rapidly as kp increases with the results summarised in figure 3. 


Rr l km! ) 
tl 13 1-5 17 ‘3 


Binding Energy / Particle (Mev) 


Fig 3: The lower curve represents the reults of a typical LOCV 
calculation including tensor correlations but excluding isobar 


effects. The upper curve results when isobars are included. 
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At last we have a curve that saturates at the empirical binding energy and 
density. An added bonus is that the curvature of the saturation curve is 
now reduced and yields a compressibility modulus K e 200-300 MeV. Further 
studies of asymmetric fluids predict a symmetry term for the semi empirical 


mass formula of 32 MeV. 


We have calculated the three-body cluster contributions to the saturation 
curve and have found them to be negligible (<1 MeV) at densities up to 
pa3po after which they grow rapidly. Thus at higher densities LOCV 
calcualtions are expected to fail and higher terms in the cluster expansion 
must be included. 


There are still many nuclear structure calculations carried out using as an 
effective interaction the BHF reaction matrix. Our analysis of the nuclear 
matter saturation problem suggests that this will underestimate the 
attraction in the *S, channe] and overestimate the attraction in the 'S, 
channel. While this latter effect is not great for pp, it is density 
dependent. It is my belief that authors who have found a need for 
effective repulsive three-body forces are simply rectifying their omission 
of the dependence of the isobar effect. 


POSTSCRIPT 


I would claim that all the essential features of the nuclear matter 
saturation problem have now been solved. The critical elements missing 
from the earlier BHF calculations were an adequate treatment of tensor 


correlations and isobar excitations of the nucleon. ¥ 


There remains, of course, a much greater mystery. With the discovery of 
the quark-gluon substructure of nucleons, isobars and mesons, why does such 


a simple potentia} mode} work at all? 


I was introduced to this problem by Hans Bethe and educated by Ben Day and 
Baird Brandow in the mysteries of the BHF approach. The LOCV calculations 
steamed from collaborations with Ray Bishop, John Owen and Majid Modarres 


with usefu] tuition on the isobar problem from Tony Green. 


Details of calculations and further reference are given in ref. [1]. 
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STRESZCZENIE 


Sewi-empiryczna formuła masowa dobrze opisuje energię wią- 
zania jąder atomowych w szerokim zakresie mas. Człon objętoś- 
ciowy w tym wzorze pozwala oszacować , przy zaniedbaniu oddzia- 
ływania Coulomta, że w jednorodnym, nieskończonym płynie skła- 
dającym się z neutronów i protonćw, tworzy się układ związany 
przy energii wiązania 16 + 5 MeV/nukleon. Ze względu na trudno- 
ści z jednoznacznym rozdzieleniem członu objętościowego i po- 
wierzchniowego, eksperymentalnie szacuje się wartość parametru 
charakteryzujacego ściśliwość płynu materii jadrowej: K =% 250 
50 MeV. Przegląd prób wyjaśnienia tej wartości jest przedmiotem 


artykułu, 


PESDME 


lonyskinupuyeckad MaccoBan dopMyna XOpomo OMMCHBaeT 3Hep- 
THO CBA3U aTOMHHX ASP B HKUDOKOM AUANA30HE Macc. OOseMinił YACH 
B 3Toli QopMyjle nosBanneT OWEHMTL, B CHiyuae NpeHGÓpeXGHKA KyO- 
HOBCKUM B3GKVMOĄCACTEMKEM, UTO B OAHOPOAHOH, OÓeCKOHEUHOŃ KMAKOCTH 
cocrazieHHOŃ u3 HeHTPOHOB M MpOTOHOB, OÓpPA3JETCA CBA3BHKAĄ CHC- 
Tesa npu 3Hepruu cBaau 1625 MaB/HykNOH, Mo nozony Tpymiocreń 
C OMKO3HAUHHM Dpa3HEJEHUEM OÓBEMHOTO X NOBEPXHOCTHOTO UHEHOB -= 
Ha OCHOBE ƏKCNEpHMEHTOB ONEHKBAETCA 3HAYEHKE IaDaMeTpa xapak- 
Tepu3ypiero cKMMaEMOCTb ZAKĄKOCTM AAEpHOŃ MaTepums K = 250£50 
iiaB. Oć3op paóoT, NATaDHKXCH BHACHATK 3TO 3HAUEHKE ABNAETCA 
npeĄueTOM 3TOŃ CTATBK. 


